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ABSTRACT 


V'  It  is  proposed  that  electromagnetic  radiation  of  a 
specified  frequency  can  be  produced  as  a  result  of  stimulated 
Cerenkov  radiation  in  a  dielectric  resonator  excited  by  a 
superluminal  electron  beam.  The  frequency  generated  is  a 
function  of  three  physical  parameters.  They  are  the 
electron  energy,  the  thickness  of  the  dielectric  resonator 
and  its  index  of  refraction.  This  work  provides  a  theoretical 
derivation  for  predicting  the  frequency  of  stimulated  Cerenkov 
radiation  in  a  dielectric  slab.  The  first  experimental 
results  using  extremely  relativistic  electrons  are  reported, 
and  the  problems  encountered  are  outlined  with  some  sugges¬ 
tions  for  improvements.  The  results  of  this  validation  show 
that  the  observed  frequency  differs  from  the  predicted  fre¬ 
quency  by  less  than  1.5%.  Incidental  to  the  conduct  of  this 
experiment,  ordinary  Cerenkov  radiation  in  the  usual  cone 
was  observed  in  air  at  microwave  frequencies.  A  possible 
application  of  the  stimulated  Cerenkov  process  as  an  electron 
beam  monitor  is  briefly  discussed.  -  - 
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I .  INTRODUCTION 


Cerenkov  radiation,  the  radiation  which  is  generated  by 
a  charged  particle  moving  at  superluminal  velocity  in  a 
medium  has  been  considered  as  a  possible  radiation  source 
for  many  years.  Additionally,  Cerenkov  radiation  has  been 
used  as  a  charged  particle  detector  which  will  selectively 
detect  only  charged  particles  exceeding  the  speed  of  light 
in  the  medium  through  which  they  are  traveling.  This  detec¬ 
tion  is  dependent  upon  only  the  velocity  of  the  particle  and 
the  index  of  refraction  of  the  medium.  The  detector  is  a 
simple  photocell  and  as  such  does  not  distinguish  between 
frequencies  within  its  band  of  operability.  Recently, 
Professor  John  E.  Walsh  (Ref.  1]  reported  successful  genera¬ 
tion  of  Cerenkov  radiation  using  relativistic  electrons  and 
a  dielectric  resonator  in  the  form  of  a  cylindrical  anulus. 
The  import  of  Walsh's  work  is  that  the  frequency  of  the 
radiation  generated  in  the  dielectric  is  a  function  only  of 
the  energy  of  the  electrons,  the  thickness  of  the  dielectric, 
and  the  index  of  refraction  of  the  dielectric.  Since  these 
quantities  are  material  parameters,  it  would  appear  that 
radiation  of  any  desired  frequency  could  be  obtained  by 
selecting  the  proper  set  of  parameters. 

The  intent  of  this  experiment  is  to  extend  Walsh's  work, 
using  a  different  geometry  for  the  dielectric  resonator  and 


■ 


much  higher  electron  energies.  Professor  Walsh  conducted  his 
experiments  at  300  KeV  in  contrast  to  the  present  experiment, 
which  uses  electron  energies  in  the  10-100  MeV  range,  an 
increase  of  1.5  to  2  orders  of  magnitude.  In  this  experiment, 
radiation  in  the  X  band  was  selected  for  ease  of  measurement 
and  analysis,  with  a  view  to  extending  the  work  into  the  IR 
range . 
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II.  THEORY 


Consider  a  dielectric  slab  with  a  coordinate  system  as 
shown  in  Figure  1.  The  dielectric  has  thickness  h  in  the  z 
direction  and  is  unbounded  in  the  x  and  y  directions.  There 
is  a  dielectric-conductor  interface  at  2=0  which  lies  in  the 
x-y  plane.  The  dielectric  slab  has  permittivity  e  and 
permeability  u o  while  these  quantities  are  e0  and  u0, 
respectively,  in  vacuum. 

Starting  with  Maxwell's  equations  and  following  the  usual 
procedure,  it  can  be  shown  that  the  longitudinal  components, 

E  or  H  ,  satisfy  a  wave  equation  and  the  modes  can  be 

X  X 

divided  into  TE  or  TM  modes.  The  general  expressions  for 
the  transverse  components  of  the  wave  equation  solution  are: 
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The  intent  of  this  experiment  is  to  take  energy  from  a  beam 
of  electrons  in  the  vacuum  just  above  the  dielectric  slab 
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and  give  this  energy  to  an  EM  wave  in  the  slab  and  in  the 
vacuum  above  the  slab.  This  energy  transfer  requires  a 
component  of  E  in  the  x-direction.  Therefore,  the  wave  in 
the  dielectric  should  be  TM  mode.  Thus,  for  TM  modes  where 


H  =0  and  E  exist,  the  wave  equation  for  E  is: 

XX  X 


(V2  +  {u2ue-k2})Ex  =  0 


(5) 


When  the  condition  that  Hx=0  is  applied  to  Equations  1,  2,  3 


V 

and  4 ,  the 

result  is 

H 
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For  simplicity,  let 

(jL)2u0£o-k~  =  -a‘ 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


for  waves  propagating  in  vacuum  and  let 

aj2UoS-k2  =  b2 

for  waves  propagating  in  the  dielectric.  Assuming  no  y 
dependence,  which  is  appropriate  for  a  slab  which  extends  to 
infinity  in  the  +  and  -  y  directions,  Eq.  5  reduces  to 


E 

— -  +  (w2ye-k2)E  =  0 

3z  2  x 


(12) 
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Thus,  in  vacuum,  substituting  Eq.  10  into  Eq.  12  yields 

3  2E 

— ^  -a2Ex  =  0  U3) 

a  Z 

which,  when  the  restriction  that  E  goes  to  zero  as  z  goes 

x 

infinity  is  applied,  has  as  its  solution 

Ex  =  Ae”az  (14) 

Similarly,  substituting  Eq.  11  into  Eq.  12  for  the  dielectric 
case  yields 

32E 

— ^  +  b2  E  =  0  (15) 

3z2  x 

When  the  boundary  condition  of  a  conductor  at  the  origin  is 
imposed,  the  solution  to  Eq.  15  is 

Ex  =  B  sin  (bz)  (16) 

At  the  vacuum-dielectric  interface,  the  usual  boundary 
conditions  apply.  That  is,  the  tangential  components  of  E 
and  H  must  be  continuous  and  the  normal  components  of  D  and 
S  must  also  be  continuous.  For  this  problem,  with  the  given 
coordinate  system,  the  tangential  components  of  £  are  E  and 
Ey,  the  tangential  components  of  H  are  Hx  and  H^.  The 
normal  component  of  6=e£  is  eE  and  the  normal  component  of 

Z* 

B=(l/u)H  is  “Hz.  Ey  and  are  zero  when  the  assumption  of 
no  y  dependence  is  applied  to  Eqs.  8  and  9  and  for  TM  modes 
Hx  is  zero  by  definition.  Thus,  the  boundary  conditions 
reduce  to:  Ex,  and  eEz  must  all  be  continuous  at  the 


vacuum-dielectric  interface.  Equating  Eqs.  14  and  16  at  z=h 


and  applying  the  boundary  condition  of  continuity  to  E  leads 

X 


Ae"ah  =  B  sin  (bh)  (17) 

Making  the  appropriate  substitutions  into  either  Eqs.  6  or  7 

and  applying  the  continuity  requirement  for  eE  at  z=»h 

z 

results  in  the  same  equation. 


icoeo(~aAe  a^) 


ioje(bB  cos{bh}) 


This  reduces  to 


£ aAe  _  eB  cos (bh) 
a  b 


Equations  17  and  19  comprise  a  coupled  set  of  equations  which 
must  be  simultaneously  satisfied. 

The  following  development  indicates  one  method  of  solu¬ 
tion  which  leads  to  an  expression  for  the  propagation 
constant,  k,  and  the  corresponding  phase  velocity  for  the 
propagating  mode.  Dividing  Eq.  17  by  Eq.  19  gives 


~  7  tan  (bh) 

to  e 


•  a  =  b  tan  (bh) 


—  ah  =  bh  tarn  (bh) 
£  o 

Let  y=ah  and  x=bh 


y  =  ~x  tan  x 
*  e 


Applying  the  definition  of  the  index  of  refraction  (n) 


y  = 


x  tan  x 
n5 


(24) 


Recalling  the  definitions  of  a2  and  b2  and  adding  them 

a2  +  b2  =  w2u0 (e-£o) 
a2h2  +  b2h2  =  u)2  u  o  ( £— £  o )  h2 
Substituting  for  (ah) 2  and  (bh) 2 
y2  +  x2  =  ui2uo  (e-£o )  h2 


(25) 

(26) 

(27) 


The  right  hand  side  of  Eq.  26  is  a  constant,  which  we  will 
call  R2 


y2  +  x2  =  R2  (28) 

The  simultaneous  solution  of  Eqs.  23  and  28  for  a  given 
frequency  will  result  in  the  value  of  k  which  yields  fields 
satisfying  all  boundary  conditions.  A  representative  solution 
of  this  system  of  equations  with  n=1.461  (polyethylene)  is 
shown  in  Figure  2. 

Although  the  TM  modes  are  the  desired  modes,  it  is 
necessary  to  also  consider  the  TE  modes  to  see  if  there  is 
any  degeneracy.  Equations  1,  2,  3  and  4  reduce  to 
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GRAPHICAL  SOLUTION 


FIGURE  2 


The  TE  mode  wave  equation  for  H  is 

X 

(72  +  {u>2iie-k2 })  H  =  0  (33) 

Again,  assuming  no  y  dependence,  Eq.  33  becomes 

32H 

37^  +  (w2ue-Jc2)Hx  =  0  (34) 

Substituting  Eq.  10  into  Eq.  34  for  waves  propagating  in 
vacuum 


3  2H 


-  a2H  =  0 


3z2  "x 

which  has  a  solution 
Hx  -  Ce-az 

and  substituting  Eq.  11  into  Eq.  34  for  waves  is  the 
dielectric 


(35) 


(36) 


3  2H 


+  b2H  =  0 


3z2 

which  has  a  solution 


H  =  D  sin  (bz) 

X 


(37) 


(38) 


when  the  appropriate  boundary  condition  at  z=0  is  imposed. 

The  same  requirements  for  the  vacuum-dielectric  interface 

apply.  The  assumption  of  no  y  dependence  immediately  makes 

Eqs.  29  and  30  zero.  The  remaining  tangential  component  of 

H  is  the  remaining  tangential  component  of  E  is  E^  and 

for  the  normal  components  of  D  and  B,  only  B  remains. 

z  v  z 

These  three  components  must  be  continuous  at  the  interface 
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at  z=h.  Thus 


Ce  =  D  sin  (bh) 


(39) 


and 

-iu>Ma  (-aCe  a^)  -iwyg  (bD  cos  { bh }) 


Dividing  Eq.  39  by  Eq.  41, 
a  =  b  tan  (bh) 
ah  =  bh  tan  (bh) 


(40) 

(41) 

(42) 

(43) 


Again,  letting  y=ah  and  x=bh, 

y  =  x  tan  x  (44) 

Equations  27  and  28  are  still  valid  and  the  simultaneous 
solution  of  Eqs .  28  and  44  will  satisfy  all  boundary  condi¬ 
tions  for  a  given  frequency. 

The  simularity  between  Eq.  44  and  Eq.  23  indicates  that 
both  TE  and  TM  modes  will  propagate  at  about  the  same  fre¬ 
quency  in  the  dielectric.  However,  only  the  TM  mode  is 
capable  of  gaining  energy  from  the  electrons.  Accordingly, 
if  stimulated  Cerenkov  radiation  is  observed,  it  will  be  the 
TM  modes  that  are  excited.  If,  however,  an  attempt  is  made 
to  amplify  waves  fed  into  the  dielectric  from  some  outside 
source,  some  mechanism  must  be  found  to  exclude  the  TE  mode 
waves  from  being  introduced  into  the  dielectric. 


A  clearer  picture  of  the  relationships  can  be  seen  by 
use  of  geometry.  First,  define 

k^  =  u2y0£o  (45) 

the  magnitude  of  the  free  space  k  vector  for  a  plane  wave 
and 

k^  =  co2Uoe  (46) 

the  magnitude  of  the  k  vector  for  a  plane  wave  in  the 
dielectric.  Substituting  these  into  Eqs.  10  and  11,  the 
defining  equations  for  a2  and  b2  respectively,  yields  after 
multiplication  by  h2, 

k2h2  -  k2h2  «  -a2h2  (47) 

k2h2  -  k2h2  =  b2h2  (48) 

Subtracting  Eq.  47  from  Eq.  48 

R2  =  h2k2  -  h2k2  (49) 

or 

R2+  h2k2  =  h2k2  (50) 

This  defines  a  right  triangle  as  shown  in  Figure  2.  Re¬ 
arranging  the  terms  in  Eq.  47  defines  another  right  triangle, 
a2h2  +  k|h2  =  k2h2  (51) 

also  shown  in  Figure  3.  Note  that  the  side  k^h  is  common 
to  both  triangles  and  that  the  magnitude  of  the  unknown  k, 
the  wave  number  of  the  propagating  mode  in  the  dielectric  that 
satisfies  all  the  boundary  conditions,  is  between  k^  and  k^. 
There  is  an  important  consequence  of  this  development. 
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Since,  for  energy  exchange,  the  phase  velocity  of  the  wave 
in  the  dielectric,  as  represented  by  the  unknown,  k,  must 
be  matched  to  the  velocity  of  the  electrons.  Also  the 
velocity  of  the  electrons  must  be  less  than  the  speed  of 
light  in  the  medium  above  the  dielectric.  Thus  the  region  above 
the  dielectric  must  be  in  vacuum  because  100  MeV  electrons 
have  a  velocity  slightly  faster  than  the  velocity  of  a 
plane  wave  in  air  under  normal  conditions. 

At  this  point  it  is  possible  to  proceed  to  a  numerical 
solution.  First,  by  phase  matching  the  velocity  of  the 
electrons  to  the  velocity  of  the  stimulated  wave  in  the 
dielectric 


Vel  vmode 


By  definition 


S  =  2 
c 


and  after  applying  Eq.  52 
~  vmode 

Again,  by  definition 


(52) 


(53) 


[54) 


U) 


mode 


mode 


mode 


U) 

3c 


(55) 

(56) 


and 


Y 2  = 


1-S2 


(57) 


B2  -/I  -  — 


(58) 
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Substituting  Eq.  58  into  Eq.  56 


1-  —  ^ 
mode  1  .  he 

"  yT> 

(59) 

Thus,  ^mode  is  determined  for  a  given  y  and 

w.  Once  u>  is 

specified,  Eqs.  45  and  46  specify  k^  and  k^. 

Eq.  10  and  substituting  Eq.  45  yields 

Rearranging 

=  kmode  '  kf 

(60) 

a  *  (kmode  -  kf 1  H 

(61) 

Similarly 

b!  -  lkd  '  k™de> 

(62) 

b  *  (kd  '  Koto'*' 

(63) 

Solving  Eq.  21  for  h 

~  *  b  tan  (bh) 

£  o 

(21) 

T7b  =  tan  (bh) 

bh  -  tan-1  (n2g) 

(64) 

(65) 

h  =  g  tan”1  (n2^) 

(66) 

Thus,  given  a  desired  frequency,  u>,  the  index  of  refraction 
of  the  dielectric,  n,  and  the  energy  of  the  electrons,  y, 
the  thickness  of  the  dielectric  can  be  determined. 

The  substitution  of  Eq.  59  into  Eqs .  10  and  11  allows  the 
factoring  out  of  w2  and  Eq.  66  becomes 
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1 


which  can  be  solved  for  ^  given  h,  n,  and  y. 


The  significance  of  Eq.  68  is  that  the  physical  parameters 
of  h,  n,  and  y  specify  oj.  With  the  proper  selection  of  these 
parameters,  EM  radiation  of  any  frequency  can  be  generated. 

Examination  of  Figure  2  readily  shows  that  if  the  dimen¬ 
sionless  parameter  bh  exceeds  rr,  more  than  one  solution  is 
possible.  However,  these  modes  will  be  at  different  frequen¬ 
cies  and  the  graphical  solution  indicates  that,  if  x<3tt,  the 
single  desired  frequency  can  be  extracted  by  use  of  a  band¬ 
pass  or  high  pass  filter. 
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III.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 


It  was  decided  to  conduct  this  experiment  in  the  X-band 
(8-12  GHz)  after  consulting  with  Professor  Knorr  and  investi¬ 
gating  the  availability  of  the  necessary  equipment  and 
measuring  devices  needed  to  support  the  experiment.  The 
dielectric  chosen  was  polyethelyne  which  has  an  index  of 
refraction  of  1.461  [Ref.  2].  The  Naval  Postgraduate  School's 
Linear  Accelerator  (LINAC)  is  capable  of  producing  relativistic 
electrons  with  energies  up  to  120  MeV.  The  operating  charac¬ 
teristics  of  the  NPS  LINAC  can  be  found  in  Appendix  A. 

A  TI-59  programmable  calculator  program  was  developed  to 
solve  Eqs .  67  and  68.  An  explanation  of  this  program  and 
a  listing  of  it  can  be  found  in  Appendix  B.  Table  I  shows  the 
input  parameter  values  and  the  program's  solution  of  Eq.  66 
for  the  0 ,  1  and  2  modes . 


TABLE  I 

THICKNESSES  OF  POLYETHELYNE  FOR  MODES  0,1,2 


m 

E  [MeV] 

n 

f  [FHz] 

h  [mm] 

0 

50 

1.461 

10 

.095 

1 

50 

1.461 

10 

14.07 

2 

50 

1.461 

10 

28.15 

^Program  output  value 


Based  on  these  results,  a  thickness  of  12.7iran  (0.5  in)  was 
selected.  The  values  of  Es50  MeV,  n=1.461,  hsl2.7mm  for  a 
Mode  1  solution  predict  a  frequency  of  11.08  GHz  when  the 
TI-59  program  solves  Eq.  68. 

The  beam  of  the  LINAC  is  approximately  1  cm  in  diameter. 
The  width  of  the  dielectric  was  made  to  be  5.08  cm  to 
approximate  the  infinite  slab  in  the  y  direction.  The  end 
of  the  dielectric  away  from  the  LINAC  exit  window  was  tapered 
to  fit  into  an  X-band  hollow  waveguide  and  to  facilitate  the 
TM  to  TE  transition  required  since  a  TE  mode  is  the  lowest 
mode  that  will  propagate  in  a  hallow  waveguide.  Figure  4  is 
a  photograph  of  the  dielectric  resonator  as  it  was  originally 
configured. 

The  first  attempt  to  observe  the  desired  EM  radiation  met 

with  only  limited  success.  The  dielectric  slab  was  mounted 

about  1  cm  below  the  centerline  axis  of  the  LINAC  exit 

window  in  air.  The  metal  waveguide  was  sloped  slightly  to 

allow  clear  passage  of  the  electrons  over  the  top  of  the 

flange  connector  to  the  detector.  A  diode  detector  was  used 

with  its  output  lead  to  an  oscilliscope  (CRO) .  With  a  beam 

2 

energy  of  53.85  MeV  and  an  average  beam  current  of  8.0  nAmp, 
the  CRO  display  showed  a  pulse  of  about  1  usee  with  strength 
varying  from  100-200  mV  as  the  beam  was  tuned,  focused  and 

2 

All  currents  read  on  a  secondary  emission  monitor  with 
6%  efficiency,  and  the  monitor  current  is  reported  in  this 
work. 
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varied  in  distance  from  the  slab.  However,  attempts  to 
measure  the  frequency  with  a  high  Q  calibrated  cavity 
resonator  were  unsuccessful. 

To  determine  if  the  observed  pulse  was  caused  by  the 
dielectric,  the  slab  was  removed  leaving  the  hollow  wave¬ 
guide  in  place.  When  the  beam  was  turned  on,  a  0.5  micro¬ 
second  pulse  of  slightly  reduced  magnitude  was  observed. 

At  this  point,  it  was  realized  that  the  observed  pulse 
might  be  caused  by  Cerenkov  radiation  in  air.  Although 
commonly  accepted  as  1.0,  the  actual  index  of  refraction 
of  standard  dry  air  is  1.003  [Ref.  3].  Using  this  value 
as  it  can  be  easily  shown  that  electrons  with  energies 

exceeding  20.856  MeV  will  exceed  the  speed  of  light  in 
standard  dry  air.  Thus,  electrons  above  that  energy  would 
exceed  the  velocity  of  a  plane  wave  in  air,  and  could  not  be 
matched  to  the  TM  mode  velocity? 

A  simple  experiment  was  designed  to  investigate  this 
possibility.  The  equipment  arrangement  is  shown  in  Figure  5. 
The  reflector  used  was  a  thin  sheet  of  aluminum  which  would 
allow  the  passage  of  the  electrons  and  reflect  any  EM  radia¬ 
tion.  With  the  reflector  in  place,  a  pulse  was  observed  from 
detector  2  whose  magnitude  was  approximately  0.5  times  the 
magnitude  of  the  output  from  detector  1.  When  the  reflector 
was  removed,  a  pulse  of  magnitude  0.1  times  the  output  of 

^See  Section  II. 
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IR  CERLNKOV  EXPERIMENTAL  SETUP 


detector  1  was  observed  from  detector  2.  These  results 
support  the  hypothesis  of  air  Cerenkov.  Since  the  air 
Cerenkov  should  be  very  broadband,  no  attempt  was  made  at 
this  time  to  measure  the  frequency  of  this  radiation. 

Based  on  these  results,  it  was  decided  to  enclose  the 
dielectric  slab  in  a  vacuum  chamber.  Figure  6  is  a  photo¬ 
graph  of  the  interior  of  the  vacuum  chamber.  The  vacuum 
was  maintained  in  the  chamber  by  installing  a  1/16  in.  thick 
piece  of  polyethelyne  between  two  waveguide  flanges  exterior 
to  the  chamber.  A  waveguide  run  of  about  75  ft.  was 
added  to  the  setup  which  allowed  for  the  positioning  of  the 
detector  in  the  LINAC  control  room  thus  eliminating  the  long, 
lossy  coaxial  cable  runs.  A  Tektronic  491  Spectrum  Analyzer 
(491-S/A)  •  •  obtained  to  facilitate  frequency  analysis. 

The  use  this  new  setup  with  team  energy  of  100.97  MeV, 
beam  current  of  4.3  nA  and  the  same  diode  detector  used  before 
resulted  in  a  0.5  Volt  peak  pulse  of  about  1  microsecond 
duration.  The  increase  in  peak  output  was  attributed  to 
the  reduction  in  line  loss  provided  by  the  replacement  of  the 
coax  cable  by  waveguide.  Subsequent  tests,  which  will  be 
discussed  later,  have  modified  this  hypothesis.  The  491-S/A, 
with  a  20  dB  attenuator  inserted  between  the  waveguide  to 
coax  adapter  and  the  S/A,  detected  signals  at  12.39  GHz  and 
8.57  GHz.  The  12.39  GHz  signal  was  the  image  of  a  real 
signal  above  the  frequency  range  of  the  491-S/A.  Using  a 
relatively  low  Q  cavity,  it  was  determined  that  most  of  the 
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power  in  the  pulse  was  in  the  8.57  GHz  signal.  A  dispersive 
high  pass  filter  with  a  cutoff  frequency  of  9.9  GHz  was 
installed  between  the  waveguide  from  the  dielectric  and  the 
coax  adapter.  With  this  arrangement  and  no  attenuation  at 
the  input  of  the  491-S/A,  valid  signals  were  detected  at 
8.58,  8.99  and  11.42  GHz.  When  the  dielectric  was  removed 
from  the  vacuum  chamber  and  the  test  repeated,  the  same 
signals  were  detected.  Both  with  and  without  dielectric,  all 
signals  disappeared  when  electron  injection  was  interrupted 
with  LINAC  RF  power  maintained. 

The  8.57  GHz  and  11.48  GHz  signals  appear  to  be  the 
third  and  fourth  harmonics  of  the  LINAC  operating  frequency. 
The  exact  cause  of  the  8.99  GHz  signal  is  not  known  but  it 
is  believed  to  be  caused  by  a  resonant  mode  of  the  vacuum 
chamber.  The  chamber  is  a  rectangular  aluminum  box  with 
10  mil  aluminum  windows  for  electron  entry  and  exit  and  as 
such  would  act  as  a  cavity  resonator.  Equipment  limitations 
prevent  the  measurement  of  any  change  in  strength  of  the 
11.48  GHz  signal  with  and  without  the  dielectric  slab  in 
place.  This  frequency  is  only  4%  above  the  predicted 
frequency  for  the  given  parameters  and  it  is  possible  that 
the  desired  signal  is  hidden  by  the  LINAC  harmonic  at  that 
frequency.  Also,  the  actual  index  of  refraction  of  the 
polyethylene  was  not  experimentally  determined  and  it  was 
not  manufactured  specifically  for  optical  use.  It  is 
possible  that  a  slight  change  in  n  and/or  inhomogeneities 
in  the  material  could  cause  some  frequency  shifting. 
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Additional  tests  were  conducted  with  this  configuration. 
Two  anomalies  were  observed  which  required  further  explana¬ 
tion.  First,  an  intermittent,  valid  signal  was  detected  at 
10.24  GHz.  This  signal  is  not  fully  understood  but  it  is 
believed  to  be  extraneous  to  the  experiment.  Secondly,  the 
peak  output  of  the  diode  detector  varied  as  much  as  an  order 
of  magnitude  with  the  same  beam  parameters  and  beam  position. 

This  observation  led  to  a  more  intensive  investigation 
of  the  TM  to  TE  coupling  at  the  dielectric/waveguide  transi¬ 
tion.  In  particular,  for  the  m=l  mode,  it  was  realized  that 
cancellation  of  the  E-fieid  could  occur  at  the  transition 
point.  A  graphical  representation  of  the  E-field  in  the 
dielectric  is  shown  in  Figure  7.  It  was  decided  to  nullify 
the  lower  field  at  the  transition  by  inserting  a  vertical 
conductor  into  the  dielectric  slab  at  the  mouth  of  the  wave¬ 
guide  as  shown  in  Figure  8.  The  vertical  conductor  should 
not  affect  any  TE  modes  propagating  in  the  dielectric  but 
should  nullify  the  cancellation  effect  of  the  m=l  mode  for 
any  TM  mode  waves  propagating  in  the  dielectric.  Another 
modification  of  the  setup  was  also  made  at  this  time.  The 
10  mil.  aluminum  window  between  the  LINAC  and  the  vacuum 
chamber  was  replaced  by  a  10  mil.  mylar  window  to  minimize 
beam  spreading  at  the  entrance  to  the  vacuum  chamber. 

Table  II  shows  the  results  of  the  frequency  analysis  of 
the  dielectric  with  and  without  the  conductor  in  place.  The 
beam  parameters  were  as  follows:  E=65.29  MeV,  1=4.0  nA. 
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VERTICAL  CONDUCTOR  IN  DIELECTRIC 


TABLE  II 


FREQUENCIES  OBSERVED 


With 

Conductor 

Without 

Conductor 

f  [GHz] 

Strength 

f  [GHz] 

Strength 

8.57 

very  strong 

8.57 

strong 

9.00 

moderate 

9.00 

weak 

9.52 

weak 

— 

— 

10.21 

moderate 

10.25 

moderate 

10.92 

weak 

— 

— 

11.42 

strong 

11.50 

moderate 

36 


The  peak  voltage  observed  from  the  diode  detector  was  150  mV 
without  the  vertical  conductor  and  100  mV  with  the  conductor 
in  place.  It  is  interesting  to  note  that  although  half  of 
the  area  of  the  waveguide  was  blocked,  the  detector  voltage 
dropped  by  only  about  one  third. 

For  input  values  of  E=65.29  MeV,  n=1.461,  h=12.7  mm  and 
m=l,  the  TI-59  program  predicts  a  frequency  of  11.08  GHz. 

The  9.52  GHz  signal  does  not  correspond  to  any  predicted 
frequency.  Possible  explanations  of  this  frequency  will  be 
discussed  later.  The  observed  10.92  GHz  signal  is  only 
1.4%  below  the  predicted  value.  Since  this  frequency  is  only 
observed  when  the  vertical  conductor  is  in  place,  it  is 
believed  that  this  is  a  TM  mode  resulting  from  stimulated 


Cerenkov  radiation. 


IV.  DISCUSSION  AND  CONCLUSIONS 


It  appears  from  the  foregoing  experiment  that  stimulated 
Cerenkov  radiation  can  be  produced  at  the  predicted  frequency 
with  the  proper  choice  of  parameters.  The  1.4%  error  can  be 
explained  in  several  ways.  First,  the  predicted  frequency 
is  based  on  an  algorithm  that  provides  an  approximate  solu¬ 
tion  for  modes  where  m>0.  Second,  nonoptical  quality  poly¬ 
ethylene  was  used  for  the  dielectric  slab.  Finally,  the 
theoretical  development  assumed  an  infinite  slab  in  the  y 
direction.  In  fact,  a  finite  slab  was  used  whose  y  dimension 
was  about  5  times  the  diameter  of  the  focused  beam.  Any  one 
of  the  deviations  from  the  ideal  case  or  a  combination  of 
them  could  account  for  the  error.  The  most  probable  source 
of  error  is  the  index  of  refraction  of  the  polyethylene.  An 
index  of  refraction  of  1.473,  a  change  of  0.012  which  is  less 
than  1%  of  the  tabulated  value  in  Ref.  3,  would  reduce  the 
error  to  essentially  zero. 

The  cause  of  the  TM  mode  frequency  at  9.52  GHz  shown  in 
Table  II  is  not  fully  understood.  It  could  be  a  result  of 
the  edge  effects  caused  by  the  finite  width  of  the  dielectric 
or  it  might  be  a  resonant  TM  mode  of  the  vacuum  chamber. 

A  secondary  effect  observed  during  the  course  of  this 
experiment  is  worth  noting.  When  the  focused  beam  was 
passed  above  the  dielectric,  the  output  of  the  diode  detector 
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varied  with  the  beam  current.  This  effect  was  present  both 
in  air  and  in  vacuum.  It  may  be  possible  to  construct  a 
calibration  curve  of  diode  voltage  vs.  beam  current  and 
thereby  develop  a  beam  current  measuring  device  which  does 
not  destroy  the  beam  downstream  from  the  measurement  point 
as  the  beam  current  measuring  device  presently  in  use  does. 

It  is  readily  apparent  from  the  magnitude  of  the  diode 
pulse  that  the  stimulated  Cerenkov  radiation  effect  is  very 
inefficient.  To  match  the  mode  with  high  speed  electrons 
requires  that  the  component  of  the  mode's  E-field  parallel 
to  the  electron  beam  be  small  compared  with  the  field's 
transverse  components.  This  will  tend  to  make  the  coupling 
small.  The  output  power  of  the  diode  detector  is  less  than 
0.5  mWatts.  Most  of  this  power,  conservatively  estimated  at 
95%,  is  in  the  8.57  GHz  signal  as  demonstrated  in  Section  III. 
This  does,  however,  enhance  the  possibility  of  developing  an 
electron  beam  current  meter,  based  on  this  process,  which 
has  only  a  minimal  effect  on  the  beam  itself. 

The  results  of  this  experiment  support  the  hypothesis 
that  stimulated  Cerenkov  radiation  can  be  produced  at  a 
frequency  specified  by  the  physical  parameters  of  electron 
energy,  index  of  refraction  and  dielectric  thickness.  It 
appears  that  the  wave-electron  coupling  is  a  weak  effect 
which  is  further  degraded  by  the  TM  to  TE  coupling  used  to 
extract  the  signal  in  this  experiment.  Greater  efficiency 
may  be  possible  if  a  more  effective  means  of  extracting  the 
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signal  is  developed.  The  signal  itself  may  be  increased  in 
strength  by  lengthening  the  dielectric  slab  to  increase  the 
length  of  the  wave-electron  interaction  region. 

Additional  work  in  the  area  of  stimulated  Cerenkov 
radiation  appears  to  be  warranted.  The  development  presented 
here  is  not  guaranteed  to  be  unique.  There  may  be  other 
modes  that  can  be  stimulated  in  the  dielectric,  one  of  which 
might  provide  a  stronger  signal.  A  mode  with  a  larger 
longitudinal  component  could  increase  the  electron-mode 
coupling.  Additionally,  the  effect  of  the  pulsed  beam  of 
electrons  should  be  studied  in  depth.  Specifically,  the 
question  of  why  this  experimental  setup  produces  most  of  its 
output  power  at  the  third  harmonic  of  the  LINAC  operating 
frequency  should  be  investigated  in  greater  detail.  If  the 
power  in  the  LINAC  harmonics  can  be  shifted  to  the  Cerenkov 
frequency  it  would  greatly  increase  the  usefulness  of  the 
stimulated  Cerenkov  process  as  a  radiation  source.  Finally, 
generation  of  X-band  radiation  using  the  m=0  mode  should  be 
attempted  as  a  preliminary  step  in  extending  this  work  into 
the  IR  range.  The  dielectric  thickness  predicted  by  the 
theoretical  development  presented  here  for  an  IR  generator 
is  of  the  order  of  microns.  The  X-band  m=0  thickness  is  in 
this  range.  Extracting  a  signal  from  a  dielectric  this  thin 
may  cause  some  difficulties.  The  experience  thus  gained 
should  be  valuable  in  extending  this  work  into  the  IR  range. 


APPENDIX  A 


LINAC  OPERATING  CHARACTERISTICS 


Max  Energy 
Max  Average  Current 
Operating  Frequency 
PRF 

Pulse  Duration 

Nr.  of  Klystrons 

Peak  Output  Power  per  Klystron 


120  MeV 
3-5  microAmp 
2.856  GHz 
60  pps 

1 . 0  microseconds 
3 

21  MWatts 
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APPENDIX  B 


TI-59  PROGRAM  LISTING  AND  EXPLANATION 

This  program  does  not  use  any  library  programs  within 

the  algorithm.  The  inputs  are  the  electron  beam  energy 

(in  MeV) ,  the  index  of  refraction  of  the  dielectric 

(dimensionless) ,  either  the  desired  frequency  (in  Hz)  or 

the  dielectric  thickness  (in  meters) ,  and  the  mode 

(dimensionless) .  The  mode  input  determines  which  leg  of 
2 

the  y=(l/n  )x  tan  x  in  the  curve  (in  the  first  quadrant) 
is  used.  For  a  given  m«0,l,2...  the  solution  is  found  for 
x  values  such  that  miT<x<  [  (2m+l)  tr/2]  .  As  noted  in  Chapter  II, 
if  x>rr,  more  than  one  solution  is  possible.  This  program 
solves  only  for  the  solution  in  the  imr<x<  [  (2m+l)  tt/2 ] 
range. 

The  solution  for  m=0  is  unique  and  exact.  Solutions  for 

higher  order  modes  are  found  using  an  iterative  search. 

The  solution  is  approximate  in  that  the  program  searches 

2 

down  the  designated  leg  of  y=(l/n  )x  tan  x  in  discrete  steps 
until  it  finds  a  y  value  less  than  the  y  value  for  the 
corresponding  m=0  case.  This  value  of  y  is  then  taken  as 
the  solution  point  for  the  input  mode. 

To  solve  Equation  67: 

1.  Load  program. 

2.  Enter  E [MeV] ,  press  A. 
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3. 

Enter 

n,  press  B. 

4. 

Enter 

f  [nz] ,  press  C. 

5. 

Enter 

m,  press  E. 

The  program  will  solve  for  h  [meters] . 

To  solve  Equation  68: 

1.  Load  program. 

2.  Enter  E [MeV] ,  press  A. 

3.  Enter  n,  press  B. 

4.  Enter  h[m] ,  press  D. 

5.  Enter  m,  press  E'. 

The  program  will  solve  for  f[Hz]. 

The  program  listing  is  on  the  following  pages. 


I 


000 

76 

LBL 

049 

65 

X 

098 

42 

STO 

001 

11 

A 

050 

89 

TT 

099 

18 

18 

002 

57 

ENG 

051 

95 

= 

100 

91 

R/S 

003 

55 

T 

052 

42 

STO 

101 

76 

LBL 

004 

93 

• 

053 

05 

05 

102 

18 

Cl 

005 

05 

5 

054 

91 

R/S 

103 

01 

1 

006 

01 

1 

055 

76 

LBL 

104 

42 

STO 

007 

01 

1 

056 

14 

D 

105 

05 

05 

008 

95 

s 

057 

42 

STO 

106 

71 

SBR 

009 

42 

STO 

058 

06 

06 

107 

88 

DMS 

010 

00 

00 

059 

91 

R/S 

108 

71 

SER 

Oil 

91 

R/S 

060 

76 

LBL 

109 

69 

OP 

012 

76 

LBL 

061 

15 

E 

110 

55 

T 

013 

12 

B- 

062 

42 

STO 

111 

02 

2 

014 

33 

x2 

063 

07 

07 

112 

55 

T 

015 

42 

STO 

064 

67 

EQ 

113 

89 

7T 

016 

01 

01 

065 

19 

D' 

114 

95 

s 

017 

65 

X 

066 

71 

SBR 

115 

91 

R/S 

018 

08 

8 

067 

88 

DMS 

116 

76 

LBL 

019 

93 

• 

068 

71 

SBR 

117 

19 

D' 

020 

08 

8 

069 

68 

NOP 

118 

71 

SBR 

021 

05 

5 

070 

43 

RCL 

119 

88 

DMS 

022 

04 

4 

071 

13 

13 

120 

71 

SBR 

023 

01 

1 

072 

65 

X 

121 

68 

NOP 

024 

01 

1 

073 

43 

RCL 

122 

91 

R/S 

025 

08 

8 

074 

11 

11 

123 

76 

LBL 

026 

52 

EE 

075 

95 

ss 

124 

10 

E' 

027 

94 

+/- 

076 

42 

STO 

125 

42 

STO 

028 

01 

1 

077 

14 

14 

126 

07 

07 

029 

02 

2 

078 

71 

SBR 

127 

67 

EQ 
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